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We analyse the role of the quark backreaction on the gauge-field dynamics and its impact on 
the Polyakov-loop potential. Based on our analysis we construct an improved Polyakov-loop 
potential that can be used in future model studies [1]. In the present work, we employe this 
improved potential in a study of a 2+\ flavour Polyakov-quark-meson model and show that the 
temperature dependence of the order parameters and thermodynamics is closer to full QCD. We 
discuss the results for QCD thermodynamics and outline briefly the dependence of our results on 
the critical temperature and the parametrisation of the Polyakov-loop potential as well as the mass 
of the (7 -meson. 
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1. Introduction 

In low-energy effective models for QCD such as the Polyakov-loop extended Nambu-Jona- 
Lasinio (PNJL) and quark-meson (PQM) model the gauge-field dynamics is not fully included. 
Specifically, the QCD Polyakov-loop potential is approximated by its Yang-Mills (YM) analogue. 
Therefore, these effective models show some discrepancy in comparison to QCD within functional 
approaches [2, 3] and lattice simulations [4-7]. To bring PNJL/PQM-type model studies closer 
to full QCD, the YM Polyakov-loop potential has to be replaced by the QCD glue potential. The 
latter is generated by the gauge degrees of freedom in the presence of dynamical quarks, for details 
see Ref. [1]. We determine the change of the Polyakov-loop potential due to the backreaction of 
the quarks on the gauge degrees of freedom. We apply this improved Polyakov-loop potential in 
a 2-1-1 flavour Polyakov-quark-meson model and compare the results with those obtained with the 
pure YM potential. A detailed discussion is given in Ref. [1], for related work see Ref. [8]. 

2. The glue potential with functional methods in pure gauge theory and full QCD 

The functional renormalisation group (FRG) equation for the scale dependent quantum effec- 
tive action, the Wetterich equation, is is depicted in Fig.l for QCD. It provides a setting which 
allows to approach the low temperature regime from large momentum scales by successively in- 
tegrating out momentum shells with momenta at about k. When lowering the momentum scale 
one systematically includes quark and gluon fluctuations into the theory finally approaching the 
hadronic phase. Although the flow equation has a simple one-loop structure. Fig. 1 describes fully 
dynamical QCD. In particular, the flow of the gluon propagator receives contributions from matter 
loops, e.g. the quark contribution to the vacuum polarisation shown in Fig. 2. 

Evaluated on a general gluonic background Aq and integrated over all scales k, Fig. 1 provides 
the full Polyakov loop potential of QCD in terms of Aq, see Ref. [2]. Here the Polyakov loop ^[Aq] 
serves as the confinement-deconfinement order parameter [9, 10]. The simple one loop structure in 
Fig. 1 allows for a straightforward identification of gluonic and matter contributions. In particular, 
the first two loops in Fig. 1 generate the glue potential of QCD and depend on the full gluon and 
ghost propagators in the presence of dynamical quarks. For example, for the gluon this entails 



Figure 1: Partially bosonised version of the FRG flow for QCD. The loops denote the gluon, ghost, quark 
and hadronic contributions, respectively. The crosses mark the RG regulator term and t = log^/A with the 
infrared momentum scale k and the references scale A. 



Figure 2: Left: Functional flow for the YM effective action. Right: Quark polarisation contribution to the 
gluon propagator representing a contribution of the matter backcoupling. 
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Figure 3: Left: YM potential (dashed lines) and glue effective potential (full lines) as functions of the 
background gauge field (Aq) for various reduced temperatures. The form of the potentials is similar, however 
the temperature scale changes. Right: Relation between the temperature scales of YM theory and the QCD 
glue potential including the backreaction of the quarks. We approximate the numerical data (black crosses) 
by a linear relation (full green line) that clearly deviates from pure YM theory (red dashed line). 

the inclusion of the quark part of the vacuum polarisation, see Fig. 2, with fully dynamical quark 
propagators. In turn, dropping the coupling between matter sector and glue sector, the first two 
loops in Fig. 1 simplify to the flow of the glue potential of pure YM theory, see Fig. 2. Within the 
FRG-approach to two flavour QCD detailed in Refs. [2, 3], we have computed the fully dynamical 
glue potential that takes into account the backreaction of the quark degrees of freedom on the 
gluon propagators (and vice versa). We compare this potential to the glue potential of SU(3) YM 
theory as obtained in Refs. [9, 11, 12]. The first non-trivial result of this comparison is that the Aq- 
dependence of both order-parameter potentials V ((Aq)) agrees on a quantitative level, see Fig. 3. 
This agreement is achieved after a non-trivial matching of the absolute temperature scales and 
the temperature slope of the potentials. Both effects on the relative temperature dependences are 
induced by the matter fluctuations altering the ghost and gluon propagators in QCD. We exploit 
these observations to estimate how the temperature-dependence of a given Polyakov loop model 
potential for pure YM theory has to be modified to improve it towards a full QCD Polyakov loop 
potential. As a measure for the distance of the potentials we use 

/ ^/(Ao) VYM((Ao))-Vglue((Ao)) . (2.1) 

JO 

In terms of reduced temperatures 

^giue = (r-r,f'=)/r,f'= and = (r - r™) /r™ (2.2) 

the minimisation of the measure (2. 1) yields the following translation of the two temperature scales 

fYM(?glue) ~ 0.57fgiue • (2.3) 

The critical temperatures of the respective potentials in the FRG computations are rj '"'^ = 203 MeV 
[2] in two-flavour QCD in the chiral limit and T™ = 276 MeV [9, 1 1, 12]. We shall discuss below 
that the absolute temperature scale still remains a free parameter (within a small interval) in our 
model. The relation (2.3) can now serve as an import input for model studies. For the detailed 
derivation and further details we refer the interested reader to Ref. [1]. 
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3. Polyakov-quark-meson model and Polyakov-loop potential 

Effective models of QCD such as the PNJL [13-16] and the PQM [17-23] model unify 
aspects of chiral symmetry and confinement. The models are defined via the free energy of QCD 
as a function of the order parameters and thermodynamic control parameters 

n(a,as,<I>,<l>;r,At/) =t/(a,as) + lt(<I>,<l>;r)+nqq(a,as,<I>,<l>;r,Ai/) . (3.1) 

The particle content of the PQM-model is constituent quarks minimally coupled to gauge fields and 
coupled to mesons via a Yukawa-type term. The contribution of the mesonic degrees of freedom 
U (cJ, as) contains several parameters that are adjusted to meson properties, see e.g. Refs. [24, 25]. 
The values of the constants that we use to fix these parameters are given in Ref. [1]. A remaining 
uncertainty in the chiral sector is the mass of the a-meson. Within our model we identify it with 
the resonance /o (500) and we discuss its impact by varying ma in the range of (400 — 600) MeV. 
The order parameters of chiral symmetry are the light or non-strange condensate a and the strange 
chiral condensate as that can be combined to the subtracted chiral condensate Ai s. 

For the Polyakov-loop potential U (<I>,0; r) a functional form is chosen that reproduces the 
temperature dependence of the Polyakov loop expectation value and the thermodynamics of pure 
YM theory as obtained in lattice calculations. We compare results obtained with the logarithmic 
parametrisation of Ref. [16] (that we abbreviate with 'Log' in the following) and the polynomial 
parametrisation with the parameters of Ref. [26] (Poly-I) and those of Ref. [14] (Poly-II). We use 
our relation between the pure gauge and the glue effective potential (2.3) to effectively convert 
these pure gauge potentials to the glue potential of full QCD 

Ugiue (<!>, ^, fglue) /T^ = Uym (<J>, ^, ?ym) /T^ ■ (3.2) 

The absolute temperature scale for the case of physical quark masses (pion masses) has not been 
computed in Ref. [2]. Therefore, we consider the glue critical temperature in the model as a free 
parameter that we vary in the range 180 MeV < Tcf'"' < 270 MeV [17], for details see Ref. [1]. 

The last term of Eq. (3.1) represents the quark sector that includes the coupling to the Polyakov- 
loop variables and the mesons. For a detailed discussion of the 2-1-1 flavour model we refer to 
Refs. [1, 18,22]. 

4. Results and discussion 

The temperature dependence of the thermodynamics and the order parameters displayed in 
Figs. 4 and 5 shows that replacing the pure YM Polyakov-loop potential by the quai^k-improved 
potential leads to a significantly smoother crossover transition at vanishing density. The pressure 
is in remarkable agreement with the result from lattice calculations given the mean field nature of 
the matter sector. The amplitude of the trace anomaly or interaction measure follows the lattice 
calculations as well. We observe that the behaviour of the pressure and interaction measure in the 
transition region is still slightly steeper. Also the evolution of the chiral order parameter shows a 
smoother evolution when applying the quark-improved Polyakov-loop potential. Including fluctu- 
ations will bring the smoothening towards the lattice result as will be discussed in a forthcoming 
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reduced temperature t reduced temperature t 



Figure 4: Normalised pressure (left) and trace anomaly or interaction measure (right) approximating the 
Polyakov-loop potential with the pure gauge potential (red dashed lines) and with the quark-improved glue 
potential (full green lines), compared to the result of lattice calculations [5, 7], see Table 1. 




reduced temperature t reduced temperature t 



Figure 5: Subtracted chiral condensate (left) and Polyakov loop (right) approximating the Polyakov-loop 
potential with the pure gauge potential (red dashed lines) and with the quark-improved glue potential (full 
green lines), compared to the result of lattice calculations [4, 6], see also Table 1. 





Uym 


^glue 


Wuppertal-Budapest [4] 


HotQCD = 12 & 8 [6] 


Tc [MeV] 


168 


158 


157 


159 & 163 



Table 1: Pseudocritical temperatures for the crossover phase transition at /i/ = 0. They are determined by the 
peaks in the temperature derivatives of the subtracted condensate Ai j. For the effective model calculations 
we use the logarithmic parametrisation of the Polyakov-loop potential with a critical temperature of 210MeV 
and a sigma meson mass of 500 MeV. 



publication. The pseudocritical temperatures summarised in Table 1 show that the result with the 
quark-improved Polyakov-loop potential agrees well with the value found in lattice calculations. 
The result for the Polyakov loop is shifted to higher values in the confined phase and to lower val- 
ues in the deconfined phase when replacing the YM Polyakov-loop potential by the effective glue 
potential, smoothening the transition. This is a consequence of applying Eq. (2.3) as can be also 
seen in Fig. 3. The remaining difference can be attributed to the missing quaik-meson dynamics 
and that the Polyakov loop potential is written in terms of the variable (<I>[Ao]) while the coupling 
to the matter sector is described as a function of [(Aq)], see the discussion in Refs. [1, 3, 27]. 
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Figure 6: Normalised pressure (left) and trace anomaly or interaction measure (right) for different parametri- 
sations of the Polyakov-loop potential and compared to the result of lattice calculations [5, 7]. Dashed lines 
stand for the case of approximating the Polyakov-loop potential with the YM potential and solid lines are for 
the case of using the improved glue potential. We use a transition temperature of the Polyakov-loop potential 
of 210MeV and a sigma meson mass of SOOMeV. 





Poly-I 


Poly-II 


Log 


rj'"' [MeV] 


180 


210 


240 


270 


Uym 


158 


169 


168 




156 


168 


183 


198 


^glue 


144 


158 


158 


Uglue 


152 


158 


164 


171 



Table 2: Pseudocritical temperatures in MeV for the crossover transition at pLf = Q. Left: For different 
parametrisations and parameter sets of the Polyakov-loop potential using a transition temperature of the 
potential of 210MeV. Right: For different critical temperatures of the Polyakov-loop potential using the 
logarithmic parametrisation. For both cases we use a sigma meson mass of SOOMeV. The pseudocritical 
temperatures are determined by the peaks of the chiral susceptibility dls.\ ^/dT . 



The uncertainties within our model study are the parametrisation and the transition tempera- 
ture of the Polyakov-loop potential and the mass of the a-meson. Figure 6 shows the sensitivity on 
the parametrisation and different parameter sets of the glue sector. The results within the different 
descriptions of the YM or glue potential are compatible but the transition is in general smoother 
with the quark-improved glue potential. The amplitude of the trace anomaly is in better agreement 
with lattice calculations. Note, however, that the temperature dependence is slightly steeper inde- 
pendent of the actual parametrisation. The maximum of the interaction measure above Tc with the 
Poly-II YM Polyakov-loop potential corresponds to the nearly constant behaviour above when 
the glue Polyakov-loop potential is applied. Table 2 summarises the pseudocritical temperatures 
for the different calculations. Both the steepness as well as the flat behaviour above Tc are expected 
shortcomings of the present mean-field nature of the matter dynamics. The effects of including a 
fully dynamical matter sector will be discussed in a forthcoming publication. 

Our analysis of the impact of the transition temperature T^"'^ of the Polyakov-loop potential 
on thermodynamic observables shown in Fig. 7 confirms that the temperature dependence of the 
pressure and the interaction measure is in far better agreement with lattice calculations when the 
quark-improved Polyakov-loop potential is used. The transition gets steeper when the ti'ansition 
temperature of the glue potential is decreased. We find best agreement with lattice calculations 
for 

T^cr ~ 2 10 MeV but this scale could be lowered by the inclusion of fluctuations [23]. The 
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8- Wuppertal-Budapest, Borsanyi etal., JHEP 11, 2010 




reduced temperature t reduced temperature t 

Figure 7: Normalised pressure (left) and trace anomaly or interaction measure (right) for different transition 
temperatures of the Polyakov-loop potential compared to the result of lattice calculations [5, 7]. Dashed 
lines are for the YM potential case and solid lines are for the improved glue potential. The logarithmic 
parametrisation of the Polyakov-loop potential and a sigma meson mass of 500 Me V is used. 




reduced temperature t reduced temperature t 

Figure 8: Same as the previous figure for different masses of the scalar sigma meson. The logarithmic 
parametrisation of the Polyakov-loop potential and a critical temperature of 210MeV is used. 



temperature dependence of the trace anomaly shows that the transition region expands below the 
pseudocritical temperature if the transition temperature Tcf'"*^ of the glue potential is decreased. The 
pseudocritical temperature is at the lower end of the transition region for larger critical temperatures 
of the Polyakov-loop potential. The absolute values of the pseudocritical temperatures are given in 
Table 2. 

The dependence of the results on the mass of the a-meson (Fig. 8) shows that the pseudo- 
critical temperature moves from the lower end of the transition region to the upper end when the 
a-meson mass is increased. Increasing has an opposite effect to increasing the transition tem- 
perature of the Polyakov-loop potential. The quantitative best agreement with the pressure of lattice 
calculations is obtained for the combinations {rricTa^^) of the a-meson mass and the Polyakov- 
loop potential transition temperature of (400MeV, ISOMeV), (500 Me V, 210MeV) and (600 Me V, 
250 MeV). The best agreement with lattice QCD data for the pseudocritical temperature is found 
for the combination = 500 MeV and TJ = 210MeV, as can be read off from Table 3. 
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m<j [MeV] 


400 


500 


600 




161 


168 


179 


^glue 


144 


158 


173 



Table 3: Pseudocritical temperatures for the crossover transition at /Z/ = for different masses of the sigma 
meson. They are determined by the peaks in the temperature derivatives of the subtracted condensate A] j. We 
use the logarithmic parametrisation of the Polyakov-loop potential with a critical temperature of 210MeV. 
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